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Abstract I present an elementary primer of black hole physics, including its gen- 
eral relativity basis, all peppered with astrophysical illustrations. Fol- 
lowing a brief review of the process stellar collapse to a black hole, I 
discuss the gravitational redshift, particle trajectories in gravitational 
fields, the Schwarzschild and Kerr solutions to Einstein's equations, or- 
bits in Schwarzschild and in Kerr geometry, and the dragging of inertial 
frames. I follow with a brief review of galactic X-ray binary systems 
with known black holes, stressing the QPO phenomenon in particular. 

I then discuss the evidence from AGN for the existence of super- 
massive black holes in galaxy nuclei, as well as evidence for such black 
holes in ordinary galaxy nuclei. I use the free motion of gas parcels 
to illustrate aspects of accretion disks around black holes, showing how 
to calculate energy efficiency and surface emissivity of disks, and the 
rate of black hole spin-up. I recall the primary methods for black hole 
mass determination, the correlation of black hole mass with the stellar 
velocity dispersion of its neighborhood, and implications for the origin 
of supermassive black holes. Finally, I consider the formation of pri- 
mordial black holes, and calculation of their mass spectrum at present 
in the case of scale invariant primordial inhomogeneities. 
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0.1 Introduction 

From the Newtonian standpoint a black hole (BH) is a compact object 
whose surface gravitational potential approaches the square of the speed 
of light. Thus in contrast with other astrophysical objects, the radius R 



2 



of a BH scales linearly with its mass M: 



R = GM/c 2 « 1.47(M/M ) km. 



(1) 



General Relativity (GR) requires the existence of BHs. Not only radius 
but several other BH properties scale simply with mass. For example 
the average mass density p = SM/AttR 3 scales as M~ 2 , 



so that massive BHs are tenuous and light BHs are dense. Classical BH 
physics lacks a special scale; thus only in specific astrophysical surround- 
ings are the physics of small and large BHs qualitatively different. 

I have chosen to discuss three separate categories of black holes: those 
that arise from stellar collapse, those that may come about from merg- 
ing and accretion processes and are termed super massive black holes 
(SMBH), and those which may be relics of the dense highly inhomoge- 
neous medium in the early universe (primordial black holes — PBH). 

No longer can it be claimed that "for astrophysical purposes a BH is 
a Newtonian point mass" ; GR has become indispensable to understand 
fine points of the observations relating BHs. Yet GR is considered by 
many opaque and hard to wield. Thus many model builders use GR for- 
mulae and results without really understanding where they come from. I 
try to alleviate the situation by providing in these lectures a little primer 
to GR together with some example calculations which are both easy to 
do and have important consequences. I also discuss some concrete issues 
in BH astrophysics, using whenever feasible GR results here obtained. 

0.2 How do black holes form from stars ? 

A good review is that of Brown, et al. [1]. A normal star is stable 
so long as nuclear burning in it provides thermal pressure to support it 
against gravity. But nuclear burning gradually converts the star core's 
hydrogen to helium, and for massive cores (M > 5M Q ) the helium is 
subsequently burned to carbon, and carbon to heavier elements until we 
reach the iron group. The core contracts as each type of fuel runs out 
because it temporarily loses pressure until the heating resulting from the 
contraction ignites the next type of fuel. Since nucleosynthesis of heavy 
nuclei from lighter ones raises the mean particle mass, it contributes to 
pressure loss. 

In the end the endothermic disintegration of iron-group nuclei, the 
tightest bound of nuclei, precipitates core collapse. Evolutionary calcu- 
lations [1] show that single stars with initial masses 20 — 3OM leave 
cores of more than 1.8M Q . These promptly collapse to "massive BHs" 



pfa 1.5 x 10 20 (M /M) 2 Kg m 
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(2) 
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(M > 1.8M Q ) with no optical display, but intense neutrino emission 
from neutronization. For initial mass 18 — 20M Q the core implosion 
(also accompanied by neutronization) engenders a shock wave which 
can blow part of the star out (type II supernova). A core of mass 
1.5M < M < 1.8M Q , momentarily neutron rich, is left; after deneu- 
tronization it collapses into a "low mass BH". Stars with initial masses 
10 — 18M© would seem to make supernovae which leave neutron star 
remnants. The fate of single stars in the initial mass range ~ 35 — 80M Q 
is uncertain. 

Binaries develop differently. According to various authors, a star with 
initial mass 20— 35M© in a binary can loose its H envelope by overflowing 
its Roche lobe to become a "naked" helium star, which upon going 
supernova leaves a low mass BH or a neutron star. Since a goodly 
fraction of stars occur in binaries, this is a way to make normal star- 
BH binaries, provided the BH formation does not disrupt the binary, 
e.g. by copious gravitational waves emission [2]. The neutron stars in 
binaries provide a two-step path for black hole binary (BHB) formation 
following Roche lobe overflow by the the companion star as it evolves 
to a giant. Accretion may take the neutron star's mass over the critical 
mass, whereupon it will implode to a BH. Of course, the residual normal 
star can itself go supernova forming a two-compact object binary. BH- 
neutron star binaries will be more common than neutron star-neutron 
star ones, but harder to see. Stars with M > 80M Q in binaries will 
make high mass BHs (as defined above). This is the origin of BHBs 
like Cyg-Xl whose ultimate fate is supernova explosion of the massive 
normal star and merger of the compact binary. 

Evidently late stellar evolution yields BHs. It is estimated that the 
galaxy contains some 3 x 10 8 stellar-mass BHs [3]. To discuss BHs we 



GR conceives gravity not as a force between masses but as a change 
in the geometry of spacetime due to neighboring energy. Consider Eu- 
clidean space which may be written in the familiar forms 



need full GR. 



0.3 



The gist of General Relativity 



ds 2 
ds 2 



dx 2 + dy 2 + dz 2 

dr 2 + r 2 {d'd 2 + sin 2 dip 2 ), 



(3) 
(4) 



or generically as 
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(5) 
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where the gij coefficients are called the metric coefficients. In principle 
there are six of them (gij counts the same as gji), but in the popular 
versions of the Euclidean line element, three vanish. If the six gij are ar- 
bitrary functions, then (5) does not necessarily describe Euclidean space, 
but some curved space in three dimensions, a Riemannian space. A clear 
example is the line element ds 2 = d"9 2 + sin 2 i? dip 2 for two-dimensional 
space. It describes the geometry of the surface of an ordinary ball of 
unit radius. This is a curved space (try flattening a world globe without 
distorting it) and it is so in two-dimensions. 

In special relativity we deal not just with space, but with flat four- 
dimensional spacetime, Minkowski spacetime, with metric 



If we replace the space part here by expression (5), we still deal the 
Minkowski spacetime. In fact 

3 



now with Greek indices like a and running over 0, 1, 2, 3 (x° represents 
the time coordinate, x % a space one), is still the spacetime of special 
relativity, provided we came by the g a p coefficients by transforming the 
coordinates from t,x,y,z in the line element (6). (Henceforth we use 
the Einstein convention: any product or single factor where a pair of up 
and down indices are equal means a sum over this index. Thus the last 
formula will be written ds 2 = g a $ dx a dx@). However, for a collection of 
ten arbitrary functions g a p (g a p counts the same as gp a ), the spacetime 
will be curved or (pseudo-)Riemannian spacetime. 

Einstein proposed that each curved spacetime represents a particu- 
lar gravitational field. Thus writing the equations of nongravitational 
physics, e.g. Maxwell's, in such a spacetime rather than in Minkowski 
spacetime (7) automatically takes care of the effects of gravity on the 
physics. This is now called the Einstein equivalence principle. Einstein 
got this revolutionary idea from the universality of free fall: all objects 
that fall freely in a given gravitational field from identical initial posi- 
tion and velocity move on identical worldlines regardless of differences 
in structure and composition. This fact is now tested experimentally 
to precision of one part in 10 13 ) (upcoming satellite experiments should 
improve the number by 3-4 orders). 

Contrast this with motion in an electromagnetic field: a positive par- 
ticle moves one way, a negative one another way, and a neutral particle 
with a dipole in yet another. What happens in gravity seems a big 



ds 2 = -c 2 dt 2 + dx 2 + dy 2 + dz 2 . 



(6) 




(7) 



a,/3=0 



Black holes: physics and astrophysics 



5 



coincidence from this perspective. Einstein elevated the universality of 
free fall in gravity to the status of a principle: the weak equivalence 
principle. He realized that the irrelevancy of structure and composition 
can be understood if particle worldlines are some universal trajectories 
in a curved geometry, so that gravity is no longer a force but a reflec- 
tion of the curvature of spacetime. Then he generalized the notion to 
all nongravitational physics — this is the Einstein equivalence principle. 
Obviously the two principles are consistent with each other. 

What are the universal trajectories ? In special relativity a particle 
not subject to forces moves inertially: if x a (r) is the trajectory {x, y, z, t} 
as functions of the proper time r, then 

d 2 x a 

Now, a theorem in Riemannian geometry tells us that locally any metric 
(7) with the correct signature can be rewritten as (6) by an appropriate 
change of coordinates. At different points we use different transforma- 
tions of coordinates, but always end up with the Lorentz metric in the 
new coordinates. So the equation (8), when written in terms of the coor- 
dinates for which the metric looks like (7), must describe the trajectory 
in the gravitational field. This is the geodesic equation (sum over (3, 7) 

d 2 x a dx? dx^ 

dr 2 ~ & dr dr [ ' 

with the functions Tg determined uniquely by (d a = d/dx a ) 



9licF% = -(&t9ni3 + dpQyn - d n9/h)- ( 10 ) 

Note that the "gravitational force" [r.h.s. of (9)] is velocity dependent. 

But is Newtonian gravitation grossly wrong ? Not at all. It is relevant 
for nonrelativistic motion, namely one with \dx l /dr\ <C c. This means 
the terms on the r.h.s. of (9) with a factor dx l /dr or two must be 
negligible. We are left with (no sum; nonrelativistically r — ► t) 

d 2 r { dt dt 

where the residual force is velocity independent. Now whenever g a p 
is time independent, e.g. field of the sun or a galaxy, and close to 
the Minkowski metric (6), we have from (10) T\ t = —\diQu- Thus (11) 
looks like the Newtonian equation for motion in a potential $ m provided 
g u = — 2$jv+ const. Comparing with the special relativistic metric (6) 
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we see that nonrelativistic motion is Newtonian if we can write 

ds 2 = -(c 2 + 2$ N )dt 2 + 2 9tidt dx i + ^ 9ijdx i dx j . (12) 

i ij 

Thus Newtonian dynamics fixes the form of gu, at least to first order 
in 3>at which is determined by Poisson's equation. We may rely on the 
above so long as |$at| <C c 2 ; otherwise the condition \dx % jdr\ <C c would 
break down. Newtonian dynamics gives no information about gu and 
gij save that they cannot be too large. 

Let us use (12) to compute the gravitational redshift, the reduction 
in the frequency of waves as they climb out of a gravitational potential 
well. Recall that we obtained (12) by assuming g a p is time independent. 
Imagine an oscillator (decaying atom, radar device) produces a wave 
train of sharp frequency v\ at a point xi. This means that N = v±Ati 
is the number of cycles of the wave in an interval Ati of the (proper) time 
ticked by a clock at rest at xi. But by (12) we have the relation Ati = 
(c 2 + 2<£jv(xi)) 1 / 2 Ai with the interval of t time spanned by the train. 
Thus the number of cycles can be written N = v\{(? + 2$jv(xi)) 1 / 2 At. 
Now the metric is not changing, so At is also the interval of t time 
spanned by the wave train anywhere, in particular at the destination X2. 
But wave cycles cannot get lost, so using his clock, the observer at X2 
must assign the wave a frequency V2 such that ^(c 2 + 2<J>at(x2)) 1 / 2 A£ = 
N. We thus have 

*1 (c 2 + 2^(x 2 )) 1 /2 ^ , AjR , , 



1/2 (C 2 + 2$ iV (xi)) 1 /2 



1 + A^/c 2 (13) 



with A&n = $jv(x2) — $jv(xi). Now the redshift z suffered by the wave 
is defined by 1 + z = A2/A1 = v\jvi\ hence Einstein's celebrated result, 

z«A$tv/c 2 . (14) 



The redshift is positive (negative) if the waves propagated uphill 
(downhill) in the potential, and its magnitude is just the change in New- 
tonian potential measured in units of c 2 . Formula (14) applies to the 
frequency of any repetitive phenomenon, e.g. the rate of arrival of pulses 
from a pulsar. Applications are many. The ratio of mass to radius of 
white dwarf stars can be deduced from the redshift of the spectral lines 
using (3). Of course one has to first correct for Doppler shift coming 
from motion along the line of sight. In the approximation which we 
work in, Doppler shift and gravitational redshift are additive. Another 
example more germane to this lecture: positronium annihilation gives 
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rise to a gamma ray line at 511 KeV. But if the annihilation occurs in 
the near environment of a gravitating body, neutron star or black hole, 
those photons will be received on Earth with a lower energy, the defect 
measuring the depth of the potential well. 

How big can the gravitational redshift get ? A neutron star of 1.5M 
has a radius of ~ 10 km giving a formal surface Newtonian potential ~ 
0.22c 2 . This already calls for nonlinear corrections. A BH, being totally 
collapsed, is more extreme. We observe from (14) that V2 — ► for any v\ 
when $at(xi) — > — c 2 /2. We are obviously pushing our formulae too far 
because they were obtained from nonrelativistic arguments, and when 
<J>jv is of order unity, motion is relativistic. Nevertheless the prediction 
that the formal Newtonian potential can only reach down to — c 2 /2 turns 
out to be correct, and so is the conclusion that the redshift is then infinite 
(see below). 

The deflection of light by gravitating bodies is another famous phe- 
nomenon. Unlike the gravitational redshift, it does depend on gu and 
gij. The reason is that photons move at speed c, so the terms in Eq. (9) 
involving T\j and T l - k are no longer small compared to that with T\ t , 
which was the only one implicated in nonrelativistic motion. T\j and 
Fj k both involve spatial derivatives of gij, so gij must be known to cal- 
culate light deflection. Newtonian arguments cannot help us here. 

We now turn to Einstein's full gravitational equation. There being 
ten metric components, there are ten partial differential equations to 
determine them. One is a fanciful elaboration of Poisson's equations 
with the relativistic energy density — as opposed to rest mass density — 
as source. Pressure and energy fluxes become the sources of the others. 
If we are mostly interested in the external gravitational field of a spher- 
ically symmetric body, then the sources can be dropped and the unique 
exact solution is Schwarzschild's metric (not Martin Schwarzschild but 
his dad Karl Schwarzschild, also the father of photographic photometry): 

ds 2 = -(c 2 + 2<Z> N )dt 2 + (l + 2c- 2 <S> N r 1 dr 2 + r 2 (d'& 2 + sm 2 tfd(p 2 ) (15) 

As usual, $tv = —GM/r where M is the mass of the object and r that 
radial coordinate for which 4nr 2 is the area of a r = const, surface. 

0.4 The black hole concept in astrophysics 

Schwarzschild's metric describes the exterior of any spherical mass- 
energy distribution as well as the simplest BH. Obviously the gu of 
this metric is the same as for the approximate metric (12). Hence the 
discussion about gravitational redshift of a source at rest goes as before, 
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and we again have 

v x 1 (c 2 + 2cI» JV (x 2 )) 1 /2 



- 1 (16) 



u 2 (c 2 + 2$ JV (x 1 ))V2 

But now we can trust the result even when is not small on scale c 2 . 
We confirm that if a source lies at r = r^, where <3?tv = — c 2 /2, then the 
redshift of radiation it emits is infinite. Thus the surface at radius 

r h = 2GM/c 2 = 2.94 M/M km, (17) 

and more generally any surface where gu = 0, is called the surface of 
infinite redshift. A fixed light source there is invisible from a distance 
because its radiation gets shifted to infinite wavelengths. 

BHs and other compact objects are frequently surrounded by accre- 
tion disks in which gas confined close to a plane moves about the object 
in nearly circular orbits. With the Schwarzschild BH in mind we can 
now understand the phenomenon of asymmetric iron line profiles in the 
X-ray spectra of BHBs. Some of these like GRS 1915+105 and V4641 
Sgr observed with BeppoSaX [4] show the unmistakable Fe Ka emission 
line (from the 2P — IS transition in Fe and its ions) with a broad red 
tail. Now, if the disk is not face on, the Doppler shift and relativistic 
beaming will convert the narrow line into two, the blueshifted compo- 
nent being the stronger. The second order (relativistic) Doppler effect 
in conjunction with the gravitational redshift will shift all frequencies 
down by a fraction 0(v 2 /c 2 ) (because v 2 /c 2 ~ <J>jv/c 2 ). But light from 
the inner parts of the disk is shifted more strongly than that from fur- 
ther out, thus also broadening the line. The three effects together thus 
broaden asymmetrically about the line's center. As we shall mention, 
the Schwarzschild accretion disk can extend down to r = 6GMc~ 2 so 
that formula (16) attests to a respectable gravitational redshift range 
quite sufficient to explain some of the mentioned observations. In other 
cases, like XTE J1650-500 observed with XMM Newton [5], a rotating 
BH (see next section), whose accretion disk can reach further in, seems 
to be called for. 

We notice that the Schwarzschild g rr metric component blows up at 
r = rh. This does not signal physically harsh conditions there, but only 
that the coordinates t and r do not retain their intuitive meaning at 
points inside r = r^. For instance, in that region the path •& = const., 
ip = const, and t = const, with r varying has ds 2 = —c 2 dr 2 < 0, 
so that it represents a physically possible motion. But who ever saw 
motion with time at a standstill ? Thus in the interior region t is not 
time, that role being usurped by r. In these lectures we do not care 
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about the interior, so we can use metric (15) with the usual meanings 
provided we are careful at r = r^. The surface r = r^ where g rr blows 
up is the boundary of the BH, and is generally called the horizon. In 
the Schwarzschild case the roles of infinite redshift surface and horizon 
are played by the same surface, but this is not a fast rule. 

The horizon can be crossed by physical entities only inward. To see 
this look at the curve i? = const., (p = const, and r = const. It has 
ds 2 = -c 2 dr 2 = -(c 2 + 2§ N )dt 2 . With r = r h , dr = 0, so this is the 
track of a "hovering" photon, and a photon is forever on the lightcone. 
But the lightcone can only be crossed inward. We conclude that in 
Schwarzschild spacetime gravity bends and blends the local light cones 
at r = rh into a spherical surface — the horizon. Anything inside it is in 
limbo for it cannot send signals out. In fact, no characteristic trace of it 
remains. A large number of calculations have made it clear that when a 
nonrotating electrically neutral object totally collapses, the BH, when it 
settles down, is a Schwarzschild BH with M as its only parameter. All 
other information about the object, its chemical composition, its thermal 
state, etc. is veiled from the exterior. A BH is the prefect shredder. 

Now we apply all this to elucidate the structure of accretion disks near 
Schwarzschild BHs. Ignoring fluid and magnetohydro dynamic effects, we 
consider free particles in circular orbits on the equatorial plane of the 
Schwarzschild metric. A key question is what is the relation between 
the azimuthal frequency v a and the radius r of the orbit. For example, 
the frequencies of the quasiperiodic oscillations (QPOs) seen in some 
galactic X-ray sources are thought to reflect the v a of close-in orbits. 

An easy procedure to v a {r) is as follows. We set r a constant and 
•d = 7r/2 (circular orbit in the plane). The r component of Eq. (9) gives 

<)%r^)%2r^)(^o. « 18 , 

The factor 2 comes about because T r tlf = Now since the metric (15) 
is (p and t independent, and g ar oc 5^, it is evident from Eq. (10) that 
T r tip = 0, while T r tt /T r w = d r g tt /d r g w . Thus u a = (27r)^ 1 Sl where 

n = <bp/dt = (-ryr^) 1 / 2 = {-d r g tt id T g w fi 2 . (19) 

Since sini? = 1, a simple calculation using Eq. (15) gives 

, , uGM\i/2 4 /M\i/lkm\§ , . 

u a = (2,)-\-^) =5.79xl0 4 (— )^(^-) 2 Hz. (20) 

Coincidentally, this has the the same form as the Newtonian azimuthal 
frequency. Other (usually smaller) frequencies associated with the orbit 
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are the radial epicyclic frequency k (the frequency of a radial perturba- 
tion of Eq. (9) about the circular orbit just discussed), and the vertical 
frequency is± (of perturbations off the exactly circular planar orbit). The 
orbit must lie at r > r^; in addition all circular orbits with < r < 3r^ 
are known to be unstable: k — > at r = 3r^ which is thus known as the 
innermost stable circular orbit (ISCO). Periodic phenomena associated 
with circular motion must take place outside the ISCO. Thus Eq. (20) 
implies for the frequency of a periodic phenomena, 

v < 2.2(M /M)kHz. (21) 

Interestingly enough, QPOs in galactic X-ray sources have frequencies 
from a few Hz reaching up to half a kHz, corresponding to the largest 
frequency allowed by (21) for a BH with 1.5 up to a few solar masses, or 
to a solar mass neutron star but with r well outside ISCO (neutron star 
radii are a few times 77J. These high frequencies are one more proof of 
the existence of compact objects in nature: even for a white dwarf v a 
for r outside the star is way too low to fit many of the observed QPOs. 

0.5 Rotating black holes 

The Schwarzschild BH does not rotate. BH physics tells us that a 
collapsed neutral rotating star gives a Kerr BH. Its line element [6] is 
parametrized by mass M and angular momentum per unit mass a: 

ds 2 = g tt dt 2 + 2g tv dtdip + g w dip 2 + EA^dr 2 + Y>d'd 2 (22) 



g tt = -(c 2 -2GMrS- 1 ) (23) 

gtip = -2GMac- 2 YT x r sin 2 d (24) 

g w = [( r 2 + a 2 c~ 2 ) 2 -aV 2 A sin 2 ^E" 1 sin 2 (25) 

S = r 2 + aV 2 cos 2 (26) 

A = r 2 -2GMc~ 2 r + a 2 C - 2 . (27) 



Note that the element gt v no longer vanishes. The Kerr parameter ac _1 
has dimensions of length. The larger the ratio of this scale to GMc~ 2 
(the spin parameter a* = ac/GM), the more aspherical the metric. 
Schwarzschild's BH is the special case of Kerr's for a = 0. 

The infinite redshift surface, this time oblate, is still where gtt = 0: 

r = roo (tf) = GMc~ 2 + [(GMc~ 2 ) 2 - a 2 c~ 2 cos 2 0] 1 / 2 . (28) 

Again, radiation reaching a distant observer from an emitter at rest there 
has its frequency gravitationally redshifted to zero. 

The horizon, that surface which cannot be crossed outward, is delin- 
eated by the condition g rr — ► oo. It lies at r = where 

r h = GMc- 2 + [{GMc- 2 ) 2 - aV 2 ] 1 / 2 . (29) 
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Indeed, the track r = r^, $ = const, with dtp/dr = a{rh 2 + a 2 ) -1 dt/dr 
has dr = (it represents a photon circling azimuthaly on the horizon, 
as opposed to hovering at it). Hence the surface r = is tangent to the 
local lightcone. The horizon meets the infinite redshift surface at the 
poles, but is otherwise inside it. Eq. (29) should not fool us; the (Boyer- 
Lindquist) coordinates used in Eq. (22) make the horizon look spherical, 
but it is not because the angular part of the Kerr metric cannot be put in 
the form f(r){dd 2 + sin 2 $ dip 2 ) reflecting the full symmetry of a sphere. 
However, all properties of the metric, including the horizon's shape, are 
axially symmetric, and of course, stationary. The horizon radius is 
well defined only for a* < 1; a BH's angular momentum has a maximum 
value that rises as M 2 . 

Inside the so called ergosphere, < r < r 00 (i?), no particle, whether 
free-falling or propelled, can stay at fixed r, <p and such track would 
have ds 2 > 0. At the very least the particle must circulate constantly 
in the angular directions. This reflects the phenomenon of "frame drag- 
ging", common to all axially symmetric metrics with gttf, ^ 0. Contrary 
to expectations from Mach's principle, the local Lorentz inertial frames 
are rotating with respect to the inertial frame at infinity (the frame of 
the stars), and when this phenomenon gets intense enough, it does not 
let the particles stay in one place. For large r, ~ r _1 while both g$$ 
and g w grow as r 2 , so that the dragging becomes rapidly imperceptible 
with growing r. It is strong only in or near the ergosphere. 

To get a better feeling for the nature of the frame dragging, let us 
find the angular velocity f2 = dip/dt of a free particle in a circular orbit 
r = const, and i9 = ir/2 about the BH. Symmetry tells us such an orbit 
will not get out of the equatorial plane. The geodesic equation (9) again 
gives Eq. (18). In view of definition (10) we now have 

{d r g tt + 2nd r9tip + Q 2 d r g w ){dt/dr) 2 = 0, (30) 

which quadratic condition for Q has the solutions 

( GM ) 1/2 (31) 

(GM)V2 oc -2 ± r 3/2 • W 



As a — > 0, the + solution here asymptotes to Schwarzschild's (19)- 
(20), as expected. The — solution obviously corresponds to a retrograde 
circular orbit since it formally gives f2 < 0. For a = both prograde and 
retrograde orbits for the same r have the same v a = (27r) _1 |f2|. But for 
o / the retrograde orbit has a larger v a than the prograde one with 
like r ! This effect, totally foreign to Newtonian theory, is a consequence 
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of the frame dragging phenomenon. With (31) goes the formula 

3.24 x W\M & /M) 
Va± = {f*T/GM)W±a. (32) 

Just as in the Schwarzschild Kerr BH has an ISCO correspond- 

ing to the vanishing of the epicyclic frequency k. (Both k and u± are 
plotted against r in Ref. [7]; they are smaller than u a .) The formulae 
for the innermost stable orbit radius nsco f° r both prograde and ret- 
rograde orbits are given by Bardeen et al. [8], and are here plotted in 
Fig. 1 together with the horizon radius r^, all as a function of a*. The 
large value of nsco ^ or retrograde orbits is what prevents the formal 
pole in Eq. (32) from showing up as a* — > 1 when — > GMc~ 2 . 



risco or r h 



0.2 0. 



Figure 1. The ISCO. The solid (broken) curve gives the radius of the ISCO of 
a Kerr BH as function of a* = ca/GM for prograde (retrograde) orbits, while the 
dotted line is the radius of the horizon. All radii are expressed in units of GMc~ 2 . 

Case + of formula (32) is plotted in Fig. 2 for the stable circular 
orbits around a Schwarzschild and extreme Kerr (a* = 1) BH. It may be 
seen that in range of r where they overlap, the curves are quite similar. 
However, for a rapidly spinning BH the curve extends much further in. 

0.6 QPOs in X-ray binary black holes 

In a detailed review of BHB, McClintock and Remillard [3] list some 40 
well observed X-ray binary sources in the galaxy with candidate BHs, 
about half of the number with high confidence. These BHBs, often 
called microquasars, exhibit five distinct spectral/temporal states, each 
source flipping among several of these: (1) the thermal-dominant (TD) 
or high/soft (HS) state, a high-intensity state dominated by thermal 
emission from the accretion disk; (2) the low/hard (LH) state, a low- 
intensity state dominated by power-law emission spectrum with a photon 
index rs 1.7, rapid variability and frequently accompanied by a radio jet; 




Figure 2. The azimuthal frequency. The azimuthal frequency v a in kHz multi- 
plied by M/Mq for prograde circular orbits outside ISCO in the equatorial plane of a 
Schwarzschild BH (dashed) and extreme (a* = 1) Kerr BH (solid curve) as function 
of rc- 2 /GM. 



(3) the quiescent state, of very faint luminosity 3 x 10 30 — 3 x 10 33 erg 
s _1 , and also dominated by power-law emission (index 1.5-2.1) whose 
very faintness is regarded as good evidence that the accreting surface is 
a horizon and no solid boundary; (4) the very high (VH) or steep power 
law (SPL) state with index 2.5 or so, and (5) the intermediate state. 
The best candidates for BHB have been seen in all of the HS, LH and 
VH states. The transition between two different states is thought to be 
related to changes in the rate of accretion to the disk. 

A respectable minority of BHB candidates exhibit bumps in the power 
spectrum of the X-ray count rate; these are called quasiperiodic oscil- 
lations (QPOs; as opposed to periodic which would imply a sharp fre- 
quency). Here we are not interested in the low frequency (0.1-few Hz) 
QPOs seen for days or weeks during the VH state with amplitudes of 
order 10% of the total X-ray counts, and whose median frequencies often 
change with the X-ray flux. Rather we shall dwell on the high frequency 
(40-450 Hz) HFQPOs seen in seven of the BHB candidates during the 
VH state. These QPOs usurp 1-3% of the X-ray flux. Three of them 
exhibit single frequency QPOs: 4U163047 (184 Hz), XTE J1859+226 
(190 Hz), XTE J1650-500 (250 Hz). Three others have each a pair of 
HFQPOs: H 1743-322 (240 Hz + 160 Hz), GRO J1655-40 (450 Hz + 
300 Hz) and XTE J1550-564 (276 Hz + 174 Hz). In addition the source 
GRS 1915+105 possesses two pairs of HFQPOs (168 Hz + 113 Hz and 
67 Hz + 41 Hz) which actually show up during the HS state. These 
QPO frequencies do not vary with the X-ray flux; they are a fingerprint 
of the system. 

The first four pairs of frequencies stand in the ratio of 3:2; in H 1743- 
322 and GRO J1655-40 the accuracy of this ratio is exquisite, in GRS 
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1915+105 it is better than 1%, and in XTE J1550-564 about 5%. Even 
more striking, the lower QPO-pair frequency has been found to follow 
the empirical law Mvq ~ 931M Hz (i/q is half that frequency and is 
identified as the "fundamental" frequency) with respect to reliable de- 
terminations of the BH masses [3] . What is this telling us ? 

A suggestive analog is resonances in our planetary system, e.g. the 
fact that the orbital period of a Saturnian satellite and that of a subset 
of the particles in the planet's rings stand in the ratio of small integers 
makes gaps in the ring. Abramowicz and colleagues [9] suggest that 
in the twin QPOs one is seeing resonances of the frequencies n and 
u± for orbits at a fixed value of rc 2 /GM for various systems. Recall 
that the Kerr metric has only one scale, GM/c 2 , and one dimensionless 
parameter, a*, and that f a 's dependence on a* is weak (Fig. 2). Orbits 
at fixed rc 2 /GM around various BHs are thus similar in a geometric 
sense, and phenomena involving just them could be expected to be alike 
despite the range of M involved. 

Now whereas in Newtonian gravity is a , k and are identical, they 
differ for orbits around (Kerr) BHs. In fact almost by definition k be- 
comes relatively small as the orbit approaches the ISCO. So, for example, 
there are stable orbits at a fixed rc 2 /GM around Kerr BHs for which 
is exactly 3/2; the resonance could amplify radial and off-plane 
perturbations, and the consequent "rattling" might well get imprinted 
on the X-ray power spectrum. Although it is not clear what the physics 
of the modulation is, the repetition of the 3/2 strongly suggests that the 
same resonance shows up in all four BHBs. Now since by Fig. 2 Mv a is 
fixed for given rc 2 /GM, and k and v± are then definite fractions of v a , 
we see that the McClintock-Remillard vq oc M _1 law is clean evidence 
that we are dealing with QPOs of BHs. 

If these ideas are correct, then the azimuthal frequency associated 
with each QPO pair must obey Mv a > 3 x O.931M kHz (recall that k 
and v± lie below the corresponding v a ). Then Fig. 2 discloses that the 
QPO source must lie at r < 5GMc~ 2 : QPO pairs are associated with 
compact objects, and they are inner disk phenomena. In fact, 5GMc~ 2 
falls below rjsco for Schwarzschild, so that some QPO pairs are Kerr BH 
phenomenon. Abramowicz et al. [9] have estimated a* ~ 0.9 for GRO 
J1655-40: rapidly rotating BHs exist in nature (more in Sec. 1.0.9). 

0.7 Supermassive black holes in galactic nuclei 

Soon after the discovery of QSOs, Salpeter [10] and Zel'dovich [11] 
suggested accretion onto a supermassive black hole (SMBH) as the QSO 
energy source. Lynden-Bell [12] solidified this understanding by stressing 
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that the large energy in a QSO radio halo, ~ 10 J, were it produced 
by nuclear reactions whose maximal efficiency is .007, would require 
the "burning" of 1O 9 M . This massive a refuse, concentrated within 
10 13 m (a dimension required by the ~day timescale of QSO X-ray 
flux variability) speaks for 10 55 J of gravitational binding energy. So 
whether nuclear power has any role in fueling a QSO or not, gravitational 
accretion power must be involved. 

The involvement of deep gravitational wells in QSO's was further ex- 
posed by the observations of seemingly "super luminal" jets emerging 
from a number of QSOs, jets which are straight for up to a Mpc. The 
"superluminar expansion" implies speeds that approach c, so the rela- 
tivists parameter GM/Rc 2 in the accelerating region must approach 
unity [13]. Rotation of a massive "gyroscope" is the logical way to sta- 
bilize the emission direction so that the jets are straight, and, of course, 
BHs can be massive rotators [14]. A good review of SMBHs in AGNs 
is given by Laor [15]. By the 1980's it was accepted that QSO's are 
ephemeral displays in nuclei of galaxies, so that if they implicate SMBHs, 
then SMBHs may be left in all galaxies which once harbored QSOs. 

More recently direct evidence for SMBHs in not obviously active 
galaxies has emerged. The evidence for some 40 of them is discussed 
in an extensive review by Kormendy and Gebhardt [16]. The strongest 
case is linked with the radio source Sgr A* at the center of the Galaxy. 
A cluster of stars is observed within 0.02 pc of it, and motions of some 
of them with velocities up to 1350 km s _1 are seen to change in time; 
the stars orbit the radio source in tens of years ! Well determined ac- 
celeration vectors for several of the stars point at Sgr A* and the data 
allow a determination of the central mass at (2.6 ± 0.2) x 10 6 M Q [17]. 
This massive object cannot itself be a cluster of smaller objects: 10 6 M 
of ordinary stars would be visible; were it brown dwarfs, they would be 
so densely packed that they would collide, merge and become luminous 
beyond observable limits. Finally a 10 6 M@ cluster of collapsed stars 
would be rapidly whittled away by evaporation [18]. 

Another strong case is in the nucleus of the spiral NGC 4258 where a 
number of H2O masers are observed. The radial (Doppler) velocities can 
be fit by rotation of a circumnuclear disk of pc scale with a Keplerian 
rotation velocity profile reaching up to 1080 km s _1 . This means there 
is a "point" mass at the center to the tune of 4 x 10 7 M Q . Once again 
we would be hard put to squeeze this massive a cluster of brown dwarfs 
or collapsed objects inside the disk's inner radius at 0.13 pc [18]. It 
must thus be a BH. The Seyfert nuclei of NGC 1068 and the Scd NGC 
4945 also have (more modest) SMBHs discovered in them by the maser 
method. 
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0.8 Supermassive black holes as energy sources 

Salpeter, Zel'dovich and Lynden-Bell all stressed the great efficiency 
with which an accreting Schwarzschild SMBH produces radiation. Later 
Bardeen [19] emphasized that BHs in nature are likely to be rotating, so 
that one should consider the energy efficiency of an accreting Kerr BH. 
Let us use the BH physics we have learned to compute the maximum 
such efficiency. The physical picture is of an accretion disk in which 
matter slowly spirals inward. As we saw, the spiraling stops at rjsco 
where an instability develops, circular motion is no longer possible, so 
that the disk ends there, if not sooner. 

First we establish that provided gas dynamics and MHD effects are 
negligible, the correct expression for relativistic total energy per unit 
mass of a gas parcel is E = —gudt/dr — gt^dtp/dr. In flat spacetime 
gu = —<? while g^p = 0, so E = c 2 dt/dr. But from special relativity 
dtjdr is the Lorentz 7, so E = 'yc 2 , indeed the relativistic energy of a 
unit mass free particle. Returning to curved spacetime we now show 
that E is conserved in free motion, so that it correctly includes the 
gravitational energy. 

In fact, whenever a metric is independent of some coordinate x^ (time 
or a spatial coordinate), the quantity K = g m dx a /dr is a constant of 
the motion for a free particle (no pressure, dissipation, etc.). For 

dK dx" dx a dx 13 dx 1 

~dr = d ? 9 ^^ ~ ^nr-dF-fr' (33) 

where we have used the geodesic equation (9) to get the second term. 
But Eq. (10) and the assumed symmetry tell us that 

9^% = ~ 2 (dp 5m7 + d i 9f3s)- (34) 

And because of the symmetry between dummy indices (3 and 7, the V 
term in Eq. (33) exactly cancels the first term. We did not need to know 
g a/ 3 explicitly to show this, only that it does not depend on x^. Since 
Kerr's metric is t independent, and gt r = gw = 0, K = —g a dt/dT — 
gtipdip/dr is conserved in free motion. By analytical mechanics it must 
be proportional to the energy, and indeed it is just E, the proposed 
expression for energy. 

In reality as a gas parcel orbits, dissipation and radiative losses cause 
its E to decrease, so that it spirals inward. Thus the E(r) for circular 
motion in Kerr as given above means the specific energy remaining by 
the time the parcel reaches r. To reduce E(r) to a practical form recall 
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that ds 2 = —c 2 dr 2 , so that for circular motion 

c 2 dr 2 = -g tt dt 2 - 2g t ^dtdv - g w d^ 2 = -g tt dt 2 (l + ^ + . 

V 9tt 9tt J 

(35) 

Substituting dt/dr from here, dcp/dt = ft as given for prograde orbits by 
Eq. (31), and the metric elements (23-27), we get after some labor 
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Figure 3. Specific energy. Ec~ per unit mass of a parcel in circular orbit of 
radius r around a Kerr BH with (in ascending order) a, = 1, 0.9, 0.8, 0. Each curve 
is shown down to the corresponding ISCO. Radii are expressed in units of GMc~ 2 . 

Fig. 3 plots Ec~ 2 for a range of radii down to rjsco for BHs with 
several spin parameters. E starts very near c 2 in the outer disk edge 
because |<E>jv| = GM/r <^ c 2 there: the parcel's energy is then its full 
rest energy. In line with our previous remark, the distribution of E 
with radius is a good predictor of how much energy is deposited at any 
particular ring of the disk; this in turn sets a bound on the emissivity. 
The steep slope of the curve for a* shows that much radiation will come 
from near the ISCO (requiring corrections for losses into the hole). 

At ISCO E reaches its lowest value, for thereafter it cannot spiral in, 
but must take a plunge into the hole so rapidly that gas dynamical pro- 
cesses make little difference, and E is very nearly conserved thereafter. 
Thus c~ 2 Ejsco signifies the fraction of the original parcel rest mass fi- 
nally accreted by the hole. Then 1 — c~ 2 Eisco is the peak efficiency 
("peak" since some of the disk's radiation may get emitted into the BH). 
This is plotted in Fig. 4 by taking Fig. 1 into account in Eq. 36. For a 
Schwarzschild BH (a* = 0) the efficiency is 0.0572, while for an extreme 
Kerr BH (a* = 1) it reaches 0.42. Thus a BH accretion disk can convert 
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Figure 4- Efficiency. The peak efficiency for conversion of mass to radiation by a 
disk around a Kerr BH with spin parameter a* . 

mass into energy up to two orders of magnitude more efficiently than 
nuclear burning: Einstein's E = mc 2 at its best ! 

0.9 The origin of black hole spin 

But why to expect a central SMBH to be surrounded by an accretion 
disk ? Gas collecting in the hole's deep gravitational well will come from 
some distance out. Whatever angular momentum it has will prevent 
it from immediately plummeting into the BH, and it will rather swirl 
around. Any rotation of the galaxy (even ellipticals and bulges rotate 
somewhat) will impose coherent swirling of different gas streams, and 
collisions among them will align the divers angular momentum vectors 
leading to formation of a thin disk. 

There is direct evidence for accretion disks in AGNs from observations 
by the X-ray satellite ASCA of the Fe Ka line coming from many such 
sources (and lately also from Chandra and XMM Newton [20]). The 
line profile shows the asymmetry discussed in Sec. 1.0.4. Models which 
assume a profile for the disk emissivity distribution consistent with the 
specific energy profile of Fig. 3 are successful in reproducing the observed 
shape of the line from the Seyfert nuclei of NGC 3516 and NGC 3227 [21]. 
In the latter case there is some indication that a rapidly rotating BH 
is involved because the strong asymmetry seen requires very relativistic 
conditions found only if the disk reaches close to the horizon. 

On this note it is easy to extend the calculations of the preceding 
section to show why, as Bardeen expected, a BH in an AGN spins rapidly. 
The physical reason is that the disk feeds the hole with a high angular 
momentum to mass ratio. To make matters simple I assume the BH is 
initially Schwarzschild with mass Mq. We need an expression for specific 
angular momentum about the axis of the disk. For nonrelativistic motion 
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in Euclidean space this would be (r sin i?) 2 dip/dt = g w £l (see Eq. (3)). 
For relativistic free motion in Kerr spacetime we should replace this by 
something similar which is conserved due to the rotational symmetry 
about the axis (as should angular momentum). The rule stated just 
ahead of Eq. (33) suggests the conserved quantity 

£ = g v tdt/dr + g vip d(p/dT. (37) 

In the limit of flat spacetime (g^t 0) this indeed reduces to the nonrel- 
ativistic expression when r — ► t, so £ is the specific angular momentum. 

When a unit mass in the disk reaches ISCO, it plummets into the 
BH without further changes to its energy E and angular momentum t. 
Thus its accretion causes changes SM = E and 5(aM) = £. Meanwhile, 
a simple calculation gives for an increment of a, = ca/GM, 

5a, = [c{GM)- 1 d{aM)/dM - 2a,]M~ 1 5M. (38) 

Putting d(aM)/dM = i/E and substituting the explicit formulae for E 
(from Sec. 1.0.8), £ from Eq. (37) and from Eq. (31), and evaluating 
at r = risco-, we obtain the differential equation for the growth of a, 
with respect to M: 

d\nM. (39) 

Here g(a,) = c 2 (GM)~ l rjsco is that complicated function of a, used to 
make Fig. 1. Integrating this equation numerically from a, = at M = 
Mq shows that the spin parameter reaches a* = 0.3058 (a* = 0.7689) 
when M has grown by 10% (40%). By the time the BH has doubled 
its initial mass, a* = 0.9875. Since a* < 1, we see that if a SMBHs 
has changed its mass substantially by accretion from its disk, it must be 
rotating rapidly. Clearly the idealization of an initially Schwarzschild BH 
is not crucial. An originally mildly rotating BH would be transformed 
into an extreme Kerr BH by the time it gained significant mass. 

0.10 How are SMBH masses measured ? 

Most SMBHs have been identified in early type galaxies [16] (excep- 
tions include the mentioned maser galaxies, the Galaxy and the two 
nearby Sbs M31 and M81). No SMBHs have been found in pure disk pr 
Ir galaxies. Detecting a BH means finding kinematic evidence (velocities 
of emission lines from a circumnuclear ionized gas disk, or velocities of 
absorption lines from bulge stars) for a large mass which is unresolved 
(more on this below). Only 10-15% of early type galaxies have well 
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formed central gas disks, but more than ten have yielded evidence for 
SMBHs, including the large E0 galaxy M87 which sports a 3 x 10 9 M 
hole as well as the famous jet which is perpendicular to the gas disk. Gas 
disk derived masses are regarded as accurate to ~ 30%. Mass determina- 
tion via stellar velocity dispersion is hampered by the long integration 
times required (ellipticals and bulges have low surface brightness), as 
well as by ambiguities of interpretation, e.g. is the velocity distribution 
anisotropic ? Yet more than 25 SMBH masses have been measured by 
this method. The range of SMBH masses is 10 6 - 10 9 M . 

The detailed tables [16] are less important than two correlations that 
emerge from them. First, the BH mass M is almost proportional to the 
bulge blue luminosity or more precisely [22, 16] 

M k 7.8 x 10 7 (L 6 /10 10 Lq) 1 - 08 . (40) 

And the BH mass is nearly proportional to the fourth power of the bulge 
velocity dispersion a e within de Vaucoleurs' effective radius [23, 16]: 

M w 1.3 x 10 8 (cr e /200km s" 1 ) 3 ' 65 . (41) 

This last has a scatter not much bigger than the mass measurement 
errors can account for (see graph in ref. [16]), so it seems to be the more 
fundamental. The SMBH accounts for about 0.0013 of the galaxy's or 
bulge's mass, but deviations of an order of magnitude are known. 

As a rule SMBHs show up in every galaxy (or bulge) observed with 
enough resolution to find a BH. In what sense ? Suppose that from the 
observed and a e an estimate is made of a central SMBH's M. Then 
the radius of its sphere of influence (region where its mass disturbs stellar 
orbits drastically) r* = GM/a^ is determined. Wherever this sphere 
can be resolved, evidence for a SMBH has always turned up. Hence 
the widely accepted belief that most early type galaxies have a central 
SMBH. SMBHs in disk galaxies with small bulges — if they exist — fail to 
comply with the M — a e relation. For example, for the Scd M33 a bound 
of 10 3 M & has been put on any central BH [16], yet the M — a e relation 
would predict one four orders heavier. SMBHs go naturally with early 
type galaxies and big bulges, not with disks. 

AGNs are generally too far for the BH sphere of influence to be resolv- 
able. For them the method of reverberation mapping has been useful. 
In a AGN there are rapid optical and UV fluctuations. Reverberation 
mapping interprets the time delays between variations in the AGN con- 
tinuum (thought to come from the accretion disk) and its broad emission 
lines (coming from enveloping gas clouds where the continuum is repro- 
cessed into lines) as the light travel times between BH and the clouds, 
and hence equal to the linear size of the region R divided by c. The 
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width of the emission lines gives a velocity dispersion a. Then, assum- 
ing that the BH is the main mass in the region one has the virial estimate 
M ~ a 2 R/G. This method provides no direct evidence for a point mass. 
And it is somewhat model dependent. Yet, after some initial confusion, 
it has been agreed that central SMBH reverberation masses lie on the 
same M — a e relation as SMBH masses in normal galaxies as determined 
by the gas disk and stellar dynamical methods. 

0.11 Where did the SMBHs come from ? 

Like questions of heritage in everyday life, the origin of the SMBHs 
abounds in controversy. Did SMBHs precede their galaxies and act as 
nuclei for their formation [24, 14] ? Or did they grow from humble 
seeds to their present proportions by accretion from an already developed 
bulge or elliptical ? Or did each develop together with its host galaxy ? 
The dearth of SMBHs in pure disks suggests that a BH is not required for 
galaxy condensation, thus militating against the first hypothesis. On the 
other hand, their presence in late-type systems is reasonable if SMBHs 
form by accretion over time; in disks rotational support would impede 
their growth, and so their absence is also reasonable. 

Early type disk galaxy bulges, like ellipticals, obey the Faber- Jackson 
relation [25] M ~ a between mass and velocity dispersion. Relation 
(41) has very much this form; hence as hinted at earlier, a SMBHs 
tends to contain a fixed mass fraction of its elliptical or bulge host. 
This again seems to support the joint growth hypothesis; for how could 
the full blown SMBH manage to control its host's final mass despite 
the small dimension of the former's sphere of influence. SMBHs have 
been found in AGNs out to z = 6, some with masses ~ 1O 9 M [26]. 
Evidently these too militate against the idea of late SMBH formation in 
fully blown galaxies. As mentioned, these early SMBHs obey the same 
rule (41) as "normal" SMBHs at the present epoch, suggesting that they 
too developed together with their hosts. 

Accordingly, Kormendy and Gebhardt [16] reach the verdict that 
AGN activity, bulge formation and SMBH gestation are facets of one 
and the same primeval process. This granted, there are remain unsolved 
problems regarding the path by which BHs grow to SMBH size, as well as 
regarding the nature of the seed BHs at the outset of the evolution. Two 
categories of growth mechanism stand out in the literature: accretion by 
BHs [27] and mergers of BHs [28]. Seed candidates include stellar-mass 
BHs, intermediate-mass (10 2 - 1O 5 M ) black holes (IMBHs) and PBHs. 
Any seed-growth method combination has to be judged by its success in 
producing the 10 9 M Q SMBHs in AGNs at z > 6. 
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A stellar-mass BH accretes slowly. Early type galaxies are gas poor, 
but have stars. To accrete a star the relatively tiny BH must first capture 
it, and then tidally disrupt it in order to ingest it, all the time respecting 
the Eddington limit. It is known that this process acting in early galaxies 
would have trouble producing 10 9 M SMBH by z = 6 [29]. But stellar- 
mass black holes could play out their role in a massive star cluster. Such 
assembly would tend to drop to the bottom of the galactic potential well 
by dynamical friction on the background, lighter, stars. Simultaneously 
distant encounters among its constituent stars would cause it to develop 
a tightly bound core and a loosely bound envelope; the last would slowly 
be lost to evaporation. 

In the dense core stellar collisions would be common and would lead 
to stellar amalgamation, and the consequent formation of massive short- 
lived stars which would evolve rapidly, some leaving BHs [30, 31]. Of 
course collisions between these last and the remaining stars should make 
the holes larger. In the meantime some of the lighter objects of either 
kind would be expelled from the cluster by Newtonian slingshot. Con- 
solidation of the cluster into an IMBH by runaway evolution would fol- 
low [31]. Other studies [32] are more pessimistic regarding the possibility 
of producing IMBHs by cluster evolution in one Hubble time. 

Elliptical galaxies (and bulges) are thought to have grown by mergers 
of smaller building blocks at z ~ 3. Being very massive on stellar scale, 
the clusters (or their IMBH remnants) would be caused to congregate at 
the merged galaxy's center by dynamical friction. Their mergers would 
thus be sped up ultimately forming baby SMBH [30]. However, this 
mechanism faces a challenge. Violent formation of a black hole can lead 
to its strong recoil by the emission of gravitational waves; this process 
can expel fresh stellar-mass BH's from the Galaxy [2]. It was recently 
realized that it can also expel newly merged SMBHs from their host 
bulges or galaxy center's [33], thus impeding their growth. 

0.12 Primordial black holes 

These problems have riveted attention on the notion that BHs origi- 
nating in the early universe are the seeds from which SMBHs grow [34] . 
PBHs, thus far conjectural, offer so much that is novel that they deserve 
a look in their own right. Because there is as yet no observational con- 
firmation of PBHs, I treat them in brief. Carr [35] and Kiefer [36] are 
recent reviews on PBH astrophysics (the second on BH physics as well). 

According to Eq. (2) the lighter a BH, the denser: a 10 9 M SMBH 
has a density like snow's; a stellar mass one is above nuclear density. 
And a BH the mass of a mountain (5 x 10 12 Kg) has a whopping density 
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2.5 x 10 55 Kgm~ 3 . Intuitively, to form a BH of such high density should 
require matter or radiation about as dense. Only in the early universe 
do we meet densities orders above nuclear density. Hence when tiny BHs 
are the issue, attention turns to cosmology, and they are termed PBHs. 

High density is not enough; a highly dense homogeneous distribution 
of matter /radiation is perfectly consistent with cosmology. To make 
BHs one needs inhomogeneities. And obviously collapse cannot take 
place until the inhomogeneity's radius R has fully entered the (particle) 
horizon (radius ct with t the cosmological time). Now during inflation 
the exponential expansion stretches every R outside the horizon (which 
grows more placidly); hence PBHs cannot form then. In addition, any 
PBHs from former eras are rapidly diluted. The post-inflation radiation 
era is the first relevant one in connection with PBHs formation. 

As an example take flat space (k = 0) radiation dominated cosmology; 
we know that the expansion factor a oc y/t. Thus the Hubble parameter 
is H = Tjt^ 1 , and the condition is R < ^cH^ 1 . But for gravitation to 
overpower pressure, R must exceed Jeans' length c s (47rGp)~ 1 / 2 , where c s 
is the sound speed and p the mass density. Due to radiation dominance 
c s = c/y/3, and since Friedmann's equation gives H 2 = (87r/3)Gp, we 
have R > ^ y /2cH~ 1 = 0.47c H -1 . Comparing we see that the collapse 
is possible only briefly after horizon crossing; thereafter R becomes very 
small compared to cH^ 1 and the second condition fails. 

From Eq. (17) for R and H « we have for the time of formation 
of a PBH of mass M (M* = 10 12 Kg): 

t » 4.9 x 1(T 24 (M/M*) s (42) 

As mentioned this is relevant for post-inflation, that is t > 10 -35 s or 
M > 1 Kg. Since the smaller the PBH, the earlier it forms, the smallest 
PBHs detectable would probe the earliest cosmological times. Once the 
PBH is formed, M is not increased sizably by accretion [35]. On the 
other hand, M can decrease by Hawking radiation. 

According to BH thermodynamics [37, 38], a Schwarzschild BH is 
endowed with a quantum temperature 

T BH = hc^iSirGkM)- 1 = 1.23 x 10 n (M*/M) °K (43) 

where k stands for Boltzmann's constant. This is manifested by the BH 
emitting into space any kind of particles in nature, each with thermal 
spectrum and statistics corresponding to the Hawking temperature Tbh- 
We see that for M 3> M*, only photons and neutrinos will be emitted. 
Anyway, this radiation is paid for by a decreasing Mc 2 . Estimating the 
energy loss with the Stefan-Boltzmann law P = (Airr^aT 4 using T 
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from Eq. (43), we get for each massless (or very light) species 

dM/dt « -4.0 x 1(T 9 (M*/M) 2 Kgs" 1 . (44) 

In reality, the BH radiation crossection is larger than 47rr 2 by a factor 
of ~ 5. And since there exist very light neutrinos and antineutrinos, 
the total dM/dt is an order of magnitude above Eq. (44). Integrating 
the fixed-up equation from M = Mq to M ~ gives the lifetime r < 
27 x 10 10 (M/M*) 3 y, the inequality appearing because as M approaches 
M*, the PBH begins emitting massive particles and M falls more quickly 
than just estimated. Evidently PBHs with M <C M* have evaporated 
by now. Nevertheless, as we shall see, PBH searches say a lot about the 
inflationary era (at whose end PBHs of only M ~ 1 Kg would form). 

PBH formation is bound up with the primordial density fluctuations 
spectrum. Inflation produces a scale invariant spectrum: the root mean 
square fluctuation e for an inhomogeneity just inside the horizon is M 
independent. When the propitious time (42) comes, a PBH can form 
only if the density contrast S = 5p/p at horizon scale is large enough; the 
rule is that for equation of state p = wp, 8 must exceed w. Evidently, 
the fraction (3 of horizon-size volumes which actually collapse (the prob- 
ability that the inhomogeneity makes a PBH), is M independent. 

To calculate the present space density of PBH formed during the 
radiation, we look at a fixed comoving volume of space; it expands as 
a 3 oc i 3 / 2 . By contrast the volume contained in a horizon is oc t 3 : the 
large volume thus contains a number of horizon-size volumes that varies 
as i -3 / 2 . Then Eq. (42) shows that the number N of PBHs of mass M 
formed out of the big volume scales like /3M~ 3 / 2 . Now the collapse takes 
about a dynamical time, GMc~ 3 , and only thereafter can the process 
repeat. Thus N refers to a range of masses ~ M so that the differential 
distribution is dN/dM oc /3M _5/2 . The power law reflects the lack of a 
special scale in the inflation spectrum. The distribution is preserved over 
time except for dilution by expansion and the dying out due to Hawking 
evaporation which only leaves PBHs with M > M*. Accordingly at the 
present epoch the differential space density is [35] 

dn/dM = l -M* 1/2 M- 5/2 Q PBH p c 9(M - M*) (45) 

where p c is the present critical density and Q PB h is the fraction of it 
in PBHs. We adjusted the numerical factor so that the integral of 
M (dn/dM) over M — the mass density — is precisely Q PBH p c . 

To calculate VL PBH focus, for example, on the end of radiation domi- 
nance, t = t r , z = 10 4 , when the horizon volume was w 4(ct r ) 2 . By to- 
day's time to that volume has expanded by a factor (10 4 ) 3 . The present 
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space density of PBHs formed then is obviously n = (5 10 _12 /(4c 3 t r 3 ). 
Eq. (45) gives the alternative estimate n rs \M J }l 2 M r ~' i l' 1 VL PBH p c . 
Equating the two, replacing M r by c 3 t r /2G, and remembering the rela- 
tion 8Gp c Ri Hq 2 between p c and Hubble constant, we get 

n pBH « 2.5 x 10 _12 /3(M r /M ;)t ) 1//2 (i7 t r ) _2 = 1.3 x 10 18 (3. (46) 

The second equality comes from two observations: for t > t r , a oc 

(2/3 

SO 

that iT = |*o and *oAr = (10 4 ) 3/2 (one computes directly M r /M* = 
6 x 10 34 ). Within the scale invariant paradigm (46) is the total PBH 
density parameter, and Eq. (45) has predictive power. 

Because we must have £I PB h < 1 we obtain a stiff inequality on (3. 
Using Carr's result (based on a Gaussian fluctuation spectrum) 

/3~ eexp(-l/18e 2 ), (47) 

we get e < 0.038. This is the accepted bound on primordial fluctuations 
from inflation [35] in the range M < 10 27 Kg. It should be stressed that 
the scale invariance assumption is a strong one, and even a small M 
dependence of e would make f3 strongly mass dependent. Thus, for ex- 
ample, Carr derives a M dependent version of Eq. 46 [Eq. (7) of Ref. [35]] , 
from which he sets constraints on e(M) separately for each mass range 
from a variety of astronomical data. Among these are the known 7 ray 
background (which would be affected by just now dying black holes), 
and the entropy in the CMB which receives a contribution from Hawk- 
ing radiation. All in all the search for PBHs makes us much wiser about 
conditions in the very early universe. 
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